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ABSTRACT 

In this paper, an analysis of an oscillatory flow of a viscous, incompressible and electrically conducting fluid with 
heat radiation in a horizontal porous channel with dissipation function is carried out. The lower stationary plate and the 
upper plate in unsteady periodic motion are subjected to a same constant injection and suction velocity respectively. 
The temperature of the upper plate in periodic motion varies periodically with time. The flow in the channel is also acted 
upon by periodic variation of the pressure gradient. A magnetic field of uniform strength is applied in the direction normal 
to the plates. A closed form solution of the problem is obtained. The effects of various flow parameters on the velocity and 
temperature fields have been shown graphically and discussed in detail. 
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1. INTRODUCTION 

The problem of finding exact solutions of the Navier-Stokes equations presents insurmountable mathematical 
difficulties. This is primarily due to the fact that Navier-Stokes equations are non-linear. This non-linearity is because of 

the presence of the convective term V.VV in these non-linear equations. There are only a few exact solutions of the 
Navier-Stokes equations known in the closed form and that too for very simple configurations of the flow patterns where 

the term V.VV vanishes in a natural way. The very basic exact solutions of the Navier-Stokes equations are found in the 
plane Couette flow, Poiseuille flow, Stokes flow etc. (Schlichting and Gersten [1]). Exact solutions for non-steady Couette 
flow were derived by Steinhauer [2] for the case when one of the walls is at rest in a steady flow and then suddenly 
accelerated to a given constant velocity. Eckert [3] obtained an exact solution of the plane Couette flow with transpiration 
cooling, which is a very effective process to protect certain structural elements in space shuttle during the re-entry, in the 
turbojet and rocket engines, exhaust nozzles etc. Singh [4] analyzed a three dimensional Couette flow with transpiration 
cooling. Taking magnetic field into account Singh and Sharma [5] have also studied MHD three dimensional Couette flow 
with transpiration cooling. A magneto hydrodynamic flow between two parallel plates with heat transfer has been analyzed 
by Attia and Kotb [6]. Chang and Lundgren [7] analyzed a duct flow in the magnetohydrodynamics. Nanda and Mohanty 
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[8] studied the hydromagnetic flow in rotating channel. Singh and Mathew [9] investigated the effects of injection/suction 
on an oscillatory hydromagnetic flow in a rotating horizontal channel. Exact solution of an oscillatory free convective 
MHD flow in a rotating channel in the presence of radiactive heat has also been studied by Singh and Garg [10]. An 
analysis of an oscillatory flow of a viscous, incompressible and electrically conducting fluid with heat radiation in a 
horizontal porous channel is studied by Singh [11]. Poonia & Chaudhary [14] and Poonia & Umashanker [15] studied 
radiation and chemical reaction effects of MHD free convective flow past an accelerated vertical plate embedded in a 
porous medium. 

The present study is aimed to analyze the oscillatory MHD flow in a horizontal porous channel with thermal 
radiation and viscous dissipation. The upper plate which is oscillating in its own plane is at a periodically varying 
temperature. 

2. MATHEMATICAL FORMULATION 


Consider the flow of an electrically conducting, viscous incompressible fluid in a horizontal channel. The two 
insulated plates of the channel are distance ‘d’ apart. The fluid is injected through the lower stationary porous plate and 
then simultaneously sucked through the upper porous plate in oscillating motion in its own plane. The constant injection 
and the suction velocities at both the respective porous plates is V. A Cartesian coordinate system (xy ') is introduced so 

that X ' -axis lies along the centerline of the channel and y ' -axis, along which a magnetic field of uniform strength B Q is 

applied, is perpendicular to the parallel plates. The magnetic Reynolds number is assumed to be very small, so that the 
induced magnetic field is negligible. The temperature difference of the plates is assumed to be high enough to induce 
radiation heat. All the physical quantities are independent of X' for this fully developed laminar flow. The flow is then 
governed by the following equations: 


cV 

a/ 


du' ,du' 

-b V - 

dt' dy' 


1 dP' d\ 

- + V- 


p du' dy' p 


°b~u r -y WU + gp T < 
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dT' ,dT' 
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1 dq' K d 2 T' 1 


dt' dy' 

dq' 
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+ - 


pC n dy' pC p dy' pC n \dy'j 


du' 


^ = 4a-r 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


(2.4) 


where, OC is the mean radiation absorption coefficient. 


The relative boundary condition can be written as 


u' = U (l + £ cos co't'),v'-V,T' = T 0 (l+ £ cos co't') at / = | 

u' = 0,v' = V,T' = 0 at y' = -~ 

2 


(2.5) 
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where, (O' is the frequency of oscillation. For the oscillatory internal flow in the channel the periodic pressure 
gradient variables are assumed to be of the form 

1 dP' 


p du' 


■ — Pcosco't' 


( 2 . 6 ) 


Because of the assumption of constant injection and suction velocity V at the lower and upper plates respectively, 
continuously equation (2.1) integrates to 


v' = V 


(2.7) 


Substituting equation (2.7) and introducing the following non-dimensional quantities 

x' y' u' 7" , , (O' P' K' 

x- — ,y = — ,u - —. 0 - — ,t - (0 t ,(0- —, P - -, K = —— 


d U T n 


v 


pUV dl 1 


into equations (2.2) and (2.3), we get 


du „ du „ dP d 2 ” ^ 


C °dt +Re dy R< dx + dy 2 
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M~+ — 
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u + GO 


( 2 . 8 ) 


n ,, 

(0 -F R, — — E 


( du^ 


dt ' dy \ dy J P dy- P 


+ - 


1 d 2 u N 


- 0 


(2.9) 


where,*, =^,P r ^ = = 

V p K K C p T {) 


sPd% 

Uv 


The boundary condition in dimensionless form become 


u — \ + £ cos t,0 — \ + £ cos t at y = - — 

2 

u-,0- 0 at y = — 

2 


( 2 . 10 ) 


3. SOLUTION OF THE PROBLEM 


For the mathematical solution of this unsteady periodic flow in the porous channel when the fluid is also acted 
upon by a periodic drop in pressure, we assume the solution in complex variable notations as 

dP 

u(y,t) = u 0 (y) + £ u i e", 6{y,t) = 0 o (y) + £0/*, ~ — = £Pe lt (3.1) 

dx 

where, P is a constant. The real part of the solution will have a physical significance. 

The boundary condition can also be written in complex notation 


www.iaset.us 


editor@iaset. us 

















52 


Umashanker, Hemant Poonia & S. S. Dhayal 


u = \ +£e“ ,6 = \ +£e u at 


u-,6- 0 at 



y = - 


2 


Substituting expressions (3.1) into equations (2.8) and (2.9), we get 


f 


U 0 R e U 0 


u\-R e u\ 
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M 1 + — 

K, 

( i A 

M 1 + — + ia> 

V K ) 


M o Gfl o 


u \ - ~PR e ~ G,B X 


0 o -P r R e 0 o -N 2 6 o =-EP r u o 

K ~ PM ~( n2 + iaP ,) 0o = -2 EPu qUi 

where, the primes in these ordinary differential equations denote differentiation w. r. t. y. 
The boundary condition (3.2) reduces to 

u 0 = 1, Mj = 1, 6 0 = 1, 6 X = 1 at y - — 

> 

u Q - 0, M; = 0 ,- 0, 0 X = 0 at y = — 

2 > 


(3.2) 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


(3.7) 


The Equations (3.3) to (3.6) are still coupled and non-linear, whose exact solution are not possible, so we can 
expand U Q ,U l ,d Q ,d l in terms of E (Eckert no.) in following form, as the Eckert number is very small for incompressible 
flows. 

u o{y)= u oi{y) +Eu o2{y) 
u i{y)= u n{y) +Eu ii(y) I 
0 o {y)=0oi{y) +E 0 O 2{y) 

0i{y)=0n{y)+ E 0i 2 {y) 


Introducing Equation (3.8) into (3.3) to (3.7), we obtain the following systems of equations. 


u oi{y) + u 'oi{y) 

K 2 {y)+K 2 {y)- 


M 1 + — 
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(3.9) 


(3.10) 


(3.11) 


Impact Factor (JCC): 4.1675 


NAAS Rating 3.45 



Unsteady MHD Poiseuille Dissipative Fluid Flow through Porous Medium with Transpiration Cooling 


u; 2 (y) + -Reu l 2 (y)- 


, 1 

M -+ — + ico 
K 


u n { y)=~ G A 


12 


0 m (y)-R e P r 0 ol (y)-N 2 0 ol (y) = 0 

- R e P Ai(y ) - N 2 0 m (y) = -P r u 0l (y ) 

0 'n(y) - PPOuiy) - (^ 2 + i(O) 0 n {y) = 0 

d n { y )~ R e P Ai [y)-(N 2 + ico)0 xl {y) = - 2 p r u 0l ( y)u n ( y) 

and the corresponding boundary conditions are 


V ^ ’ ^01 1’^02 0,65,1 1? 0 q 2 0>#n 1’ #12 ^ 

y ^01 ^02 ^12 #01 #02 #n ^ ^12 ^ 


Solving equations (3.9) to (3.16) under the boundary conditions (3.17), we get 
# 0 i (y) = a 1 e m ' > ’ + a 2 e ,niy 


u 0l (y) = a b e mj,y +a 1 e"' iy + a 8 e m ' y + a 9 e" l2 - y 
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(3.13) 
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Skin-Friction 

The expression for the shear stress is given by 

f dF A 


T — — 


\ 


d y Jy =0 


: —\u n +£u,e" ) 

V u 1 >aty=0 


Nusselt Number 

The expression for the rate of heat transfer is given by 


Nu = — 


d yjy =0 


■ ~(0 o +£0y) 

V u 1 t at y=0 


(3.26) 


(3.27) 


4. RESULTS AND DISCUSSION 

In order to have a physical insight into the problem we have evaluated, numerically the expressions for the 
velocity u(y, t ) , the Temperature 0(y,t ) , the skin friction T and the rate of heat transfer Nu. These numerical values 
are shown graphically to assess the effects of the variables of Reynolds number Re, Hartmann number M, The Pressure 
gradient P, The Porosity Parameter K, The Eckert number E and the frequency of oscillation (0. The effect of various 
parameters on skin friction and Nusselt number is also shown in table-1 and 2. 



Figure 1, shows the effects of Prandtl number Pr and the Pressure gradient P on velocity profile against y. It can 
be seen that velocity profile is increasing exponentially as increasing y. It also observed that as the favorable pressure 
gradient P increasing the velocity also increases and remains positive over the entire width of the channel. This means that 
the increasing favorable pressure gradient accelerates the flow field. It is also examined that the velocity for Pr=0.72 is 
higher than that of Pr=7.0. Physically it is possible because fluids with higher Prandtl number have high viscosity and 
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hence move slowly. 

The variation of the velocity profile with the Grashof number Gr. The frequency of oscillations CO and the Eckert 
number E are shown in figure 2. The magnitude of velocity leads to an increase with an increase in Gr. It is due to the fact 
an increase in the value of the thermal Grashof number has the tendency to increase the thermal buoyancy effect. It is also 
observed that increasing the frequency of oscillation leads to decrease in the velocity field, whereas, an increase in the 
Eckert number results an increase in the velocity field. 




Figure 3: Effects of M and K on Velocity Profile Figure 4: Effects of Re and N on Velocity Profile 
(P=5,Pr=0.72,Re=0.5,N=l,Gr=l, CO =5,E=1, £ =0.2) (P=5,Pr=0.72, M=2, K=l, Gr=l, CO =5, E=l, £ =0.2) 


Figure 3, shows the variation of velocity profiles under the influence of magnetic parameter M, the porosity 
parameter K, it is evident from figure 3 that the velocity decreases with increase of magnetic parameter. This is because of 
the reason that the effect of a transverse magnetic field on an electrically conducting fluid gives rise to a resistive type 
force (called Lorentz Force) similar to drag force and upon increasing the values of M increases the drag force which has a 
tendency to slow down the motion of the fluid. It is also examined that the velocity increases with increasing the porosity 
parameter. The presence of a porous medium increases the resistance to flow resulting in a decrease in the flow velocity. 
This behavior is depicted by the decrease in the velocity as K decreases. 

The variation of velocity profile with the Reynold number Re and the radiation parameter N are shown in figure 4. 
It is concluded that an increase in Reynold number results a decrease in the velocity profile. 

Figure 5 shows the variation of temperature under the influence of the Prandtl number Pr, the Reyonld number Re 
and the radiation parameter N against y. It is examined that temperature profile increases with increasing the Prandtl 
number for y < -0.2, after that the magnitude of temperature for air (Pr=0.72) is greater than of water (Pr=7.0). This is due 
to fact that the thermal conductivity of fluid decreases with increasing Pr, resulting a decrease in thermal boundary layer 
thickness. It also concluded that the temperature profile decreases with increasing Re. 

In both figures 4 and 5, velocity and temperature profile increase with the increase of radiation parameter N. 
The effect of radiation is to increase the rate of energy transport to the gas, there by making the thermal boundary layer 
become thicker and fluid become warmer, this enhances the effect of thermal buoyancy of the driving body force due to the 
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mass density variation which are coupled to temperature and there for, increasing the fluid velocity. 

The variation of temperature profile in the frequency of oscillations CO and the Eckert number E are exhibited in 
figure 6. It is found that the temperature field goes on increasing with the increase of CO four y < 0 and after that the 
temperature field decreases with increasing CO. It is also examined that an increase in the Eckert number leads to an 
increase in temperature profile. 




Figure 5: Effects of Pr, Re and N on Temperature Figure 6: Effects of CO and E on Temperature 
Profile, (P=5, M=2, K=l, Gr=l, CO =5, E=l, £ =0.2) Profile, (P=5,Pr=0.72, 

Re=0.5,N=l,M=2,K=l,Gr=l, £ =0.2) 


Table 1 and Table 2 shows that the effect of various flow parameters on skin-friction and Nusselt number 
respectively. It is concluded from table-1 that skin-friction increases with increasing Re, N and M, where as, it decreases 
with increasing P, Pr, K, Gr and E. It can be seen from table-2 that Nusselt number increases with increasing Re and N, 
where as, it decreases with increasing Pr and E. 

NOMENCLATURE 

U - the constant velocity of vertical porous channel ( m.S 1 ) 

T 0 - the temperature at the channel (K), T - the temperature in the boundary (K) 

U ,V - denotes the component of velocity in the boundary layer in X and y respectively ( m.S ) 
t' - the time ( S ), d - constant, g - the acceleration due to the gravity ( m.S 1 ) 

C p - the heat capacity of the fluid ( J.Kg .K 1 ), B t] - the magnetic induction, q - the Radiactive heat flux, 

P - the Prandtl number, K’ - the Porosity parameter, K - Non-dimensional Porosity parameter 
M - the magnetic parameter. Re - the Reynold number, E - the Eckert number, P - the pressure gradient 
N - the radiation parameter, G r - the thermal Grashof number, Nu - the Nusselt number 
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GREEK LETTERS 

OC - the mean radiation absorption coefficient 

(5 - the volumetric coefficient of thermal expansion ( K 1 ), p - the density of the fluid, Kg.m 

6 - Dimensionless temperature, fj, - the coefficient of viscosity (Pa.s) 

2 -1 

T - Dimensionless skin-friction, V - the kinematics viscosity (m .S ) 

<7 - the electrical conductivity, 1C - the coefficient of thermal conductivity ( Will 1 K 1 ) 

CO - the frequency of the suction velocity 


Table 1 


La. 

P 

Pr 

Re 

N 

M 

K 

Gr 

E 

-0.98132 

-10 

0.72 

0.5 

1 

2 

1 

1 

1 

-0.98985 

0 

0.72 

0.5 

1 

2 

1 

1 

1 

-0.99838 

10 

0.72 

0.5 

1 

2 

1 

1 

1 

-0.99412 

5 

0.72 

0.5 

1 

2 

1 

1 

1 

-1.17644 

5 

7 

0.5 

1 

2 

1 

1 

1 

-0.96154 

5 

0.72 

1 

1 

2 

1 

1 

1 

-0.89339 

5 

0.72 

2 

1 

2 

1 

1 

1 

-0.98783 

5 

0.72 

0.5 

2 

2 

1 

1 

1 

-0.96239 

5 

0.72 

0.5 

5 

2 

1 

1 

1 

-1.17572 

5 

0.72 

0.5 

1 

0 

1 

1 

1 

-0.63611 

5 

0.72 

0.5 

1 

4 

1 

1 

1 

-1.02814 

5 

0.72 

0.5 

1 

2 

5 

1 

1 

-1.03252 

5 

0.72 

0.5 

1 

2 

10 

1 

1 

-1.18569 

5 

0.72 

0.5 

1 

2 

1 

5 

1 

-1.3862 

5 

0.72 

0.5 

1 

2 

1 

10 

1 

-1.01755 

5 

0.72 

0.5 

1 

2 

1 

1 

5 

-1.04683 

5 

0.72 

0.5 

1 

2 

1 

1 

10 


Table 2 


Nu 

Pr 

Re 

N 

E 

-1.29427 

0.72 

0.5 

1 

1 

-8.07233 

7 

0.5 

1 

1 

-1.28016 

0.72 

1 

1 

1 

-1.20662 

0.72 

2 

1 

1 

-1.12064 

0.72 

0.5 

2 

1 

-0.25769 

0.72 

0.5 

5 

1 

-1.93679 

0.72 

0.5 

1 

5 

-2.73993 

0.72 

0.5 

1 

10 


5. CONCLUSIONS 

In this paper, we analyzed an oscillatory flow of a viscous, incompressible and electrically conducting fluid with 
heat radiation in a horizontal porous channel with dissipation function. The lower stationary plate and the upper plate in 
unsteady periodic motion are subjected to a same constant injection and suction velocity respectively. 

The velocity for Pr=0.72 is higher than that of Pr=7.0 due to fluids with higher Prandtl number have high 
viscosity and hence move slowly. It is observed that increasing the frequency of oscillation leads to decrease in the velocity 
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field, whereas, an increase in the Eckert number results in an increase in the velocity field and also examined that the 
velocity increases with increasing the porosity parameter because the presence of a porous medium increases the resistance 
to flow resulting in a decrease in the flow velocity. 

6. REFERENCES 

1. Attia, H. A. and Kotb, N. A.,(1996): MHD flow between tow parallel plates with heat transfer, Acta Mechanica, 
Vol. 117, pp-215-220. 

2. Chang, C. C. and Lundgren, T.S.,(1961): Duct flow in magnetohydrodynamics, ZAMP, Vol.l2,pp-100-l 14. 

3. Cogley, A.C.L., Vincent, W.G. and Giles, E.S.(1968): Differential Approximation for Radiating Transfer in a Non 
linear equations gray gas near Equilibrium, American Institute of Aeronautics and Astronautics, Vol.6,pp-551- 
553. 

4. Eckert, E. R. G.,(1958): Heat and Mass Transfer, McGrow Hill, New York. 

5. Nanda, R. S.and Mohanty, H. K.,(1971): Hydro magnetic flow in a rotating channel, Appl. Sci. Re s., Vol. 24,pp-65- 
78. 

6. Schlichting, H. and Gersten, K.,(2001): Boundary Layer Theory, McGraw Hill, New York. 

7. Singh, K.D.,(1999): Three dimensional Couette flow with transpiration cooling, ZAMP Vol.50, pp-661-668. 

8. Singh, K.D.and Rakesh Kumar,(2001): MHD three dimensional Couette flow with transpiration cooling, 
ZAMM, Vol. 81, pp-715-720. 

9. Job, VICTOR M., and S. Rao Gunakala. "Unsteady MHD Free Convection Couette Flow between Two Vertical 
Permeable Plates in the Presence of Thermal Radiation Using Galerkin ’s Finite Element Method." International 
Journal of Mechanical Engineering 2.5 (2013): 99-110. 

10. Singh, K.D. and Mathew, A.,(2008): Injection/suction effect on an oscillatory hydromagnetic flow in a rotating 
horizontal porous channel, Indian J. Phys. Vol.82 (4),pp-435-445. 

11. Singh, K. D. and Garg, B. P. (2010): Exact solution of an oscillatory free convective MHD flow in a rotating 
porous channel with radiative heat, Proc. Nat, Acad. Sci. India, Vol. 80(A), pp-81-89. 

12. Rao, B. Madhusudhana, G. Viswanatha Reddy, and M. C. Raju. "Unsteady MHD mixed convection of a viscous 
double diffusive fluid over a vertical plate in porous medium with chemical reaction, Thermal radiation and joule 
heating." International Journal of Applied Mathematics & Statistical Sciences (IJAMSS) 2.5 (2013): 93-116. 

13. Singh K.D. (2011): Exact solution of an unsteady periodic MHD Poiseuille flow with transpiration cooling and 
thermal radiation, Int. J. of Physical and Mathematical Sciences (IJPAMS), Vol. 2(1), pp-125-132. 

14. Poonia, H. and Choudhary, R. C. (2009): Influence of dissipative fluid on MHD free convective heat transfer flow 
through porous medium, J. of Rajasthan Acad, of Physical Sciences, Vol. 8, pp-475-484. 

15. Tiwari, A., Tiwari, K.K., Chauhan, T.S., and Chauhan, I.S., (2013): Free Convective Heat And Mass Transfer 
Flow Under The Effect Of Sinusoidal Suction With Time Dependent Permeability, International Journal Of 
Emerging Technology And Advanced Engineering, Vol.3(4),Pp-442-451. 


Impact Factor (JCC): 4.1675 


NAAS Rating 3.45 


Unsteady MHD Poiseuille Dissipative Fluid Flow through Porous Medium with Transpiration Cooling 


59 


16. Poonia, Hemant and Chaudhary, R. C. (2016): Mass transfer with chemical reaction effects on MHD free 
convective flow past an accelerated vertical plate embedded in a porous medium, Int. J. of Applied Mathematics 
& Statistical Sciences, Vol. 5, pp-33-46. 

17. Poonia, Hemant and Umashanker (2016): Radiation effect on natural convection flow past an impulsively started 
infinite vertical plate through porous medium in the presence of magnetic field and first order chemical reaction, 
Int. J. of Applied Mathematics & Statistical Sciences, Vol. 5, pp-17-28. 


www.iaset.us 


editor@iaset. us 




